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Abstract 

Motivated by the recent works of us [1, 2], we estabhsh a correspondence be- 
tween the /(i?)-gravity with the polytropic and Chaplygin gas dark energy models. 
We reconstruct the different /(i?)-gravity models according to the polytropic, stan- 
dard Chaplygin, generalized Chaplygin, modified Chaplygin and modified variable 
Chaplygin gas dark energy models. We obtain the necessary conditions for crossing 
the phantom divide line of the selected models which describe accelerated expansion 
of the universe. 
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1 Friedmann equations in the /(i?) -gravity 

The action of general /(i?)-gravity is [3] 

S^J V=^d'^{ ^^2F^^ + ^matter}, (l) 

where k'^ — SnG. Also G, g, R and ^matter are the gravitational constant, the determinant 
of the metric Qij^i,, the Ricci scalar and the lagrangian density of the matter inside the 
universe, respectively. 

Now if we consider the spatially-flat FRW metric for the universe as 

ds'' ^-dt^ + a\t)j2idx'f, (2) 
1=1 

and taking T^^ = dia.g{—p,p,p,p) for the energy-momentum tensor of the matter, then 
the Friedmann equations in the f{R) theory can be written as [4] 

^H'^ ^ Pm + PR, (3) 



^{2H + 3H')^-{p^+Pn), (4) 



where 



-^/(i?) + 3{H + H^)f\R) - 18{AH^H + HH)f"{R) 



(5) 

(6) 



\f\R)-{H + ?>H^)f\R)+^ 

6{8H'^H + 6HH + AH'^ + H)f"{R) + ?,&{HH + AH'^Hf f" {R)\ , 
with 

R^6{H + 2H^). (7) 

Here and p„i are the energy density and pressure of the matter, respectively. Also pr 
and Pr are due to the contribution of /(i?)-gravity. 
The energy conservation laws are still given by 

Pm + ^H{p^ + pm) = 0, (8) 

Pr + 3H{pr + pr)^0. (9) 
The equation of state (EoS) parameter of the curvature contribution is defined as [5] 

PR 

(^R= 

PR 

A(Hf{R) + 3(3M - AH^H + AH^ + H)f"{R) + 18{H + AHHff"{R)) 

(^-f{R) + 6{H + H^)f'{R)-36{AHm + HH)f"{R)) '^^^^ 

For a given a — a{t), by the help of Eqs. (5) and (6) one can reconstruct the /(i?)-gravity 
according to any DE model given by the equation of state (EoS) p — p{p). 
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Here we assume a pole-like phantom scale factor as [6] 

a{t) = ao{ts - ty^, t<ts, h>0. 



(11) 



Using Eqs. (7) and (11) one can obtain 

a — ao 



R 



6(2/1+1) 

h 



R, 



6h{2h+l) 

Substituting Eq. (12) into (5) and (6) gives 

1 



(12) 



PR 



2/c2(2/i+l) 



{2h + l)f{R) + {h + l)Rf\R) - 2R'f{R) 



(13) 



PR 



3h{2h+l)f{R)-h{3h+l)Rf'{R)+2{2h+3)R^f"{R)+4f"'{R) . (14) 



3k^h(2h+l) 

Also the EoS parameter of the curvature (10) yields 



Ur = -1 



2R{hf\R) -{h- 3)Rf"(R) + 2R^f"'{R) 
3h({2h + l)f{R) -{h+ l)Rf'{R) + 2R^f"{R)) ' 



(15) 



As we mentioned before to solve the above differential equations we have to know the 
EoS p — p{p) or p = p{a). In the next sections, we reconstruct different /(i?) -gravities 
according to the polytropic and different versions of Chaplygin gas dark energy (DE) 
models. 



2 Polytropic /(i^)-gravity model 

Here we reconstruct the /(i?)-gravity from the polytropic gas DE model. Following [7], 
the EoS of the polytropic gas is given by 

PA = Kpl+", (16) 

where K is a positive constant and n is the polytropic index. Using Eq. (9) the energy 
density evolves as 

PA=(5a^-i^)"", (17) 

where 5 is a positive integration constant [7]. Inserting the second relation of Eq. (12) 
into (17) one can get 

PA=iBao^{6h{2h+l))~R^-K] , (18) 
Equating (13) with (18) gives 

2R'f"{R) -{h + l)Rf{R) - {2h + l)f{R) + {e + ^R^) = 0, (19) 
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where 



-1 

„2\ " 



e = -K(2{2h+l)k^ 

— ..-in. — 1 

C = 5ao" (6/i(2/i + 1)) (2{2h + 1)A;2) " . 
Solving Eq. (19) yields the polytropic /(i?)-gravity as 
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(20) 



fp{R) = A+i?™+ + X-R" 



£"(2/i + l)Vl - lOh + h^ 



X 



, X r-, I 2nm- 2?7.m_ —^R^n 

(m+) 2-Fi I — ^7— 



3/i 



3/i ' 



(m_) 2i^i I 1 



3/i 



3/i ' 



where 2-^1 denotes the first hypergeometric function and 



m± 



(22) 



One notes that to have a real /(-R), the parameter h must be in the ranges of < /i < 0.1 
and h > 9.9. 

Replacing Eq. (21) into (15) one can get the following EoS parameter 



UJR^ -1 + 



_ „ -3h ) 

1+ 



h>0. 



Using Eqs. (12) and (20), the above relation can be rewritten as 
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UJR = -1- 



h>0. 



K 
B 



(23) 



(24) 



We see that for ^ 



do {-^1 > I5 < — 1 which corresponds to a phantom accelerating 

universe. The result of Eq. (24) is same as that obtained for the EoS parameter of the 
polytropic /(T)-gravity model given by [2]. 



3 Standard Chaplygin /(i^)-gravity model 

The EoS of the standard Chaplygin gas (SCG) DE is given by [8] 

A 

PA = , 

Pa 



(25) 



where ^4 is a positive constant. Inserting the above EoS into the energy conservation 
equation, leads to a density evolving as [8] 



Pa 



(26) 



where B is an integration constant. 
Inserting Eq. (12) into (26) gives 

PA = ^A^^EB^h, (27) 

where 

E = Bao-%6h{2h + l)f''. (28) 
Equating (27) with (13) one can obtain 



2R'f"{R) -{h + l)Rf{R) - {2h + l)f{R) + + X = 0, (29) 



Jl3h 

where 

e^A(2(2h+l)k^y, 
^ = Bao~^i6h{2h + l)f''(2{2h + 1)A;2)^ 
Solving the differential equation (29) yields 

fscciR) = A+i?'"+ + A_i?"^- X 

■' ^ ' ^ (2/i+l)Vl-10/i+/i2 



(30) 



.(31) 



Replacing Eq. (31) into (15) one can get the EoS parameter of the standard Chaplygin 
/(i?)-gravity model as 

, L... - (32) 



1 + IR^h' 



Using Eq. (30) it can be rewritten as 



B(6h(2h+1)Y'' 
= + ^ ; 73^, (33) 

which for 5 < and R > 6h{2h + l)ao^(^)^ then ujji can cross the phantom-divide 
line. 

4 Generalized Chaplygin /(i?)-gravity model 

The EoS of the generahzed Chaplygin gas (GCG) DE model is given by [9] 

PA = -4' (34) 
Pa 

where a is a constant in the range < a < 1 (the SCG corresponds to the case a — 1) 
and A a positive constant. Using Eq. (9), the GCG energy density evolves as [9] 

P'^={^ + ^) ■ (35) 
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where B is an integration constant. Inserting Eq. (12) into (35) one can get 

Pa = 



-3/j(l+a) -1 -rr— 

A + E R^^Y+", (36) 



where 

3fa(l+a) 

£; = Sao-'(^+")(6/i(2/i+l)) ' . (37) 
Equating (36) with (13), one can obtain 

2R'f{R) -{h+ l)Rf'{R) - {2h + l)f{R) + \e + ^rrl ^ = 0, (38) 



where 
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e= (2{2h + l)eyA, 

e = (2(2/i + l)k^y Bao-''\Qh{2h + 1))^, (39) 



^ — 1 + a. 

Solving the differential equation (38) yields 



fGCGiR) = A+i?™+ + A.i^"- - , = X 

^ ^ ^ + (2/i+l)Vl-10/i + /i2 

/2m_ -1 2m_ A . . ^ /2m+ -1 2m+ ' 

Replacing Eq. (40) into (15) one can get 



.(40) 



Using (39) one can rewrite the above equation as 

3fe(l+a) 

D(Gh(21i + l)) ' 

= + iMiT^, (42) 

ao3(i+")^ R^^ + B(&h{2h + 1)^ ^ 



which clears that for B < and R > 6h{2h + l)ao^ (^) '"^'^"^ then oor can cross the 
phantom-divide line. 

5 Modified Chaplygin /(i?)-gravity model 

The EoS of the modified Chaplygin gas (MCG) DE model is given by [10] 

PA = ApA-—, (43) 
Pa 
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where A and B are positive constants and < a < 1. Using Eq. (9), the MCG energy 
density evolves as [10] 

/ B C 
PA = + ^^(ir^KTM) j ' (44) 

where C is an integration constant. Inserting Eq. (12) into (44) one can get 



Pa 



B „ _ -3h(l + a)(l+A) 

+ E R 5 



1 + A 



l+a 



where 



E = Cao-'('+")('+^) {6h{2h + 1)) 
Equating (45) with (13) yields 



3h(l+Q)(l+A) 
2 



(45) 



(46) 



2R'f"{R) -{h + l)Rf{R) - {2h + l)f{R) + 



£ + 



„ 3/177) 
it! 2 



0, 



where 



£ = B(l + yl)-i (2(2/1 + 1)A;2 
^ = Cao '^''(6/i(2/i + 1))^ (2(2/i + l)k^)\ 
7 = 1 + 0;, 
77 = 1 + A 



Solving Eq. (47) gives 
/mcg(-R) = A+i?"*+ + X-R^~ — 



2e'- 



(2/i+l)Vl-10/i + /i2 



X 



(47) 



(48) 



. ^ ( 2m_ —1 2m_ — \ . s ^ f 2m+ —1 2m+ — ^ 
(m+)2Fi ( ^7^, — ; 1 + ^7:7:::; —MHE - \T^-hFi ^7 — , — ; 1 + ^7 — ; —nil 



V 3/1777 7 3/i7?7 
Replacing Eq. (49) into (15) gets 



V 3/1777 7 3/1777 sR-^ J _ 



.(49) 



= -1 + 



1 + A 



i + fit:^' 



(50) 



Using Eq. (48) it reduces to 



C(l + A)2(6/i(2/i+l)) 



aft(l + a)(l + A) 
2 



Q^3(l+a)(l+A)5 ^Ml±i^ ^ ^ j^g^^2/i + 1)) 



3>t(l + a)(l + A) ) 
2 



(51) 



2 

which shows that for C < and R > 6/z-(2/i + l)ao~ 
the phantom-divide line. 



B 



3/»(l+a)(l+A) 



then ur can cross 
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6 Modified Variable Chaplygin /(i?)-gravity model 



The EoS of the modified variable Chaplygin gas (MVCG) DE model is given by [11] 

B 



Pa = ApA - 



(52) 



where A, B and n are positive constants and < a < 1. Using Eq. (9), the MCG energy 
density evolves as [10] 



Pa = 



3(1 + 



+ 



C 



^TTJT V3(l + ")(1 +^) Q3{l+a)(l+A)-n 



where C is an integration constant. Inserting Eq. (12) into (53) one can get 



Pa = L ii:2(i+a) 



M 



+ N i?t("-3(l+a)(l+A)) 



l+a 



where 



T — n 

1j — 



M 



6/i(2/i+ 1) 

3(r+a)-B 



nh 



2(l+a) 



3(l+a)(l+A)-n' 



Equating (54) with (13) gives 



2R'f{R) -{h + l)Rf\R) - {2h + l)f{R) + 



e + 



C 1- 
RS 



= 0, 



where 



3(l+a).B[6fc(2fe+l)] ^ [2(2ft+l)fc2]^ 

[3(l+a)(l+A)-nJ ' 

C[6fe(2fe+1)]'^'^"^°''^^^' [2(2fe+l)fc2 
3(l+a)(l+A) 



7 = 1 + a, 

^ _ 3h{l+a)(l+A) _ nh 



a — 



nh 



2 ' 



Solving the differential equation (56) yields 
/mvcg(-R) — A+i2'"+ + A_i?"^" — 



2e' 



{h{2 + a) + (1 + a)(l + 2a))^/T^^WlTh? 
Replacing Eq. (58) into (15) yields the EoS parameter as 



X 



-1 + 



3/i 



l{i + eR^). 



(53) 



(54) 



(55) 



(56) 



(57) 



(58) 



(59) 
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With the help of Eq. (57) it yields 



= -1 + 



3(1 + a) 



n + 



452(7(00"' 6/i(2/^+l))^ 



(60) 



which shows that for C < and R > 6h{2h + l)ao 



'l\C\5(l+A){2nh-S) 
3nh(l+a)B 



then the EoS 



parameter of the /MvcG(-R)-gravity model corresponds to a phantom accelerating universe. 



7 Conclusions 

Here we considered the polytropic gas, SCG, GCG, MCG and MVCG models of the DE. 
We reconstructed the different theories of modified gravity based on the f{B) action in the 
spatially-flat FRW universe and according to the selected DE models. We also obtained 
the EoS parameter of the polytropic gas, standard Chaplygin, generalized Chaplygin, 
modified Chaplygin and modified varaibel /(i?)-gravity scenarios. We showed that cross- 
ing the phantom-divide line can occur when the constant parameters of the models to be 
chosen properly. 
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